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Fermat's principle of least time in the presence of uniformly 
moving boundaries and media 
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Sts. Cyril and Methodius University, P. 0. Box 162, 1000 Skopje, Macedonia 
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The refraction of a light ray by a homogeneous, isotropic and non- dispersive transparent material 
half-space in uniform rectilinear motion is investigated theoretically. The approach is an amalgama- 
tion of the original Fermat's principle and the fact that an isotropic optical medium at rest becomes 
optically anisotropic in a frame where the medium is moving at a constant velocity. Two cases of 
motion are considered: a) the material half-space is moving parallel to the interface; b) the material 
half-space is moving perpendicular to the interface. In each case, a detailed analysis of the obtained 
refraction formula is provided, and in the latter case, an intriguing backward refraction of light is 
noticed and thoroughly discussed. The results confirm the validity of Fermat's principle when the 
optical media and the boundaries between them are moving at relativistic speeds. 

PACS numbers: 03.30.+p, 42.15.-i, 02.30.Xx 



I. INTRODUCTION 

The variational approach in geometrical optics was es- 
tablished in 1657 when the famous French mathemati- 
cian Pierre de Fermat provided a mathematical proof 
that the straight line is not the fastest way for light 
to traverse between two optical media pi 0, 01 • This 
simple and new way of looking at the propagation of 
light was already empirically noticed by Hero of Alexan- 
dria fifteen centuries before Fermat by observing that 
the path of the light between two points upon reflec- 
tion from a plane mirror is the shortest possible path. 
Fermat's insightful and intriguing idea that "the Nature 
always acts in the shortest ways" has influenced gen- 
erations of scientists, like Maupertius, Euler, Lagrange, 
Hamilton and Feynman, to successively translate classi- 
cal and quantum physics into the language of variational 
calculus 0, 0, S Q ■ Today, the practical value of Fer- 
mat's principle is immense, ranging from applications in 
geophysics, oceanography and material physics on one 
side^ to general relativity and astrophysics on the other 

0, S, I M EL El m EI El El, E3, El El, IH . 

The modern formulation of Fermat's principle states 
that a light in going between two fixed points traverses an 
optical path length that is stationary with respect to the 
variations of that path. In most cases, this formulation 
reduces to the original one that the actual path between 
the two points is the path which requires the shortest 
time [2l| . An accurate description of the process beyond 
Fermat's principle is provided by the modern quantum 
electrodynamics, where the light is considered to be made 
out of photons. Accordingly, the photon takes all the 
available paths in going between two fixed points, each 
path at a certain probability, and the principle of least 



'Electronic address: 
t Electronic address: 



agj urcin@iunona. pmf . ukim.edu. mk 



skepalek@iunona.pmf.ukim.edu.mk 



time follows as an approximation of this argument [22j, 

MM- 

Fermat's principle has been initially established for 
optical systems whose properties are time-independent, 
and, as such, has entered the mainstream of the standard 
undergraduate physics curriculum through the deriva- 
tions of the well-known laws of reflection and refraction 
of light at a smooth stationary interface between two op- 
tical materials at rest [22[ . Recently, Fermat's principle 
was proven valid for optical systems with time-dependent 
parameters, incorporating the situations of li ght propa- 
gation in uniformly moving optical media [25L l26l. \2l\. It 
was also shown that Einstein's law of reflection of light 
from a plane mirror in uniform rectilinear motion fol- 
lows from Fermat's principle and the constant speed of 
light postulate, providing an indirect proof of the validity 
of Fermat's principle when the boundaries upon which 
the light is reflected are moving at relativistic speeds 
[2^ |. Apart from the latter example, it seems that the 
problems of light propagation involving uniformly mov- 
ing boundaries and media have never been attacked di- 
rectly by using Fermat's principle of least time. 

Motivated by the latest extensions to non-stationary 
cases, in this paper we will use Fermat's principle to de- 
rive the formulas for the law of refraction of a light ray in 
two specific situations when the boundaries and the me- 
dia are moving at constant velocities. We will limit our 
analysis to optical media that are homogeneous, isotropic 
and nonconducting in their rest frames of reference. Also, 
we consider the media to be non-dispersive at rest, in the 
sense that the correcting terms due to dispersion are so 
small that can be neglected. 

We begin our discussion in Sec. II by showing that 
an optical medium, isotropic in its rest frame, exhibits 
an optical anisotropy with respect to the frame where 
the medium is in uniform rectilinear motion. We derive 
the expression for the speed of the photon in the moving 
medium and show that it depends on the angle between 
the light ray and the velocity of the medium. We de- 
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scribe this relativistically induced optical anisotropy in 
a more convenient way, by introducing an effective re- 
fractive index of the moving medium, defined as a ratio 
between the speed of the photon in vacuum and the speed 
of the photon in the moving medium. The notion that 
the photon will "see" the moving medium as a station- 
ary medium having an effective refractive index different 
from its refractive index at rest, is incorporated into Fer- 
mat's principle in Sec. Ill and IV to investigate the re- 
fraction of the photon incident from a vacuum half-space 
upon a uniformly moving material half-space. In Sec. Ill 
we investigate the situation when the material half-space 
is moving parallel to the interface, and in Sec. IV we 
take the material half-space to move perpendicularly to 
the interface. In both cases, we consider the plane of 
incidence to be normal to the vacuum-material interface 
and parallel to the velocity of the medium. 

A summary of the results and a comparison with other 
similar treatments are given in Sec. V, where the possi- 
bility to extend the method to more complicated config- 
urations is also suggested. Alternative derivation of the 
refraction law when the material half-space is moving 
perpendicular to the interface is presented in Appendix 
A. Finally, in Appendix B, the conditions for the appear- 
ance of the backward refraction are derived analytically. 



II. AN OPTICAL ANISOTROPY OF A 
TRANSPARENT MEDIUM INDUCED BY ITS 
UNIFORM MOTION 



wavefront 




FIG. 1: An illustration of the light drag in the S-frame where 
the medium is moving at a constant speed u (< c/n) in the 
direction of the positive a>axis. The motion of the medium is 
a cause for the wavefront of the pulse originating from O to be 
a dragged ellipse possessing an axial symmetry with respect 
to the i-axis. Here, the photon emitted from O at an angle 
ip will travel at a speed v(tp). 



If a flash of light is emitted from a fixed point source 
in an empty space (a vacuum), the outgoing photons will 
move radially on a spherical wavefront expanding at a 
speed of light c 29] . With respect to an observer from a 
frame that moves in a straight line at a constant velocity, 
the wavefront of the pulse will remain sphere, expanding 
at the same constant speed c in all directions 30]. The 
result is a corollary of the second postulate of special 
relativity that the speed of light in a vacuum is the same 
in all inertial frames of reference regardless of the motion 
of the light source (3l| . 

The situation is changed if the expansion of the pulse 
is taking place in a material medium (32l . HH . If we con- 
sider a homogeneous, isotropic and non-dispersive trans- 
parent medium at rest with respect to 5"-frame, and if we 
treat the problem two-dimcnsionally, the equation that 
describes the evolution of the pulse is: 



+ 2 / 2 = (c/n)Y 2 , 



(1) 



where n is the refractive index of the material at rest. 
This is an equation of an expanding circle whose radius 
at time t' is (c/n)t'. By using the Lorentz transformation: 

x' — j(x — ut), y' — y, tl = j(t — ux/c 2 ), (2) 

with 7 = (1 — u 2 /c 2 ) -1 / 2 , we can make a transition to 
5-frame where the medium is moving at a constant speed 



u in the positive direction of the x-axis. In this way, Eq. 
((T|) is transformed into: 



- yt) 2 
(at) 2 (bty 

where by rj, a and b we denoted: 



I - 1/n 2 



= 1, 



(c/n) 



b = (c/n) 



1 — u 2 /n 2 c 2 / 
f - u 2 /c 2 




(3) 

(4) 
(5) 
(6) 



According to Eq. ([3]), the wavefront of the pulse (the 
ray surface) with respect to an observer in .S-frame is 
an expanding ellipse "dragged" by the moving medium 
(see Fig. 1). At time t from the emission of the pulse, 
the product rjt denotes the dragging parameter, that is, 
the distance from the origin O of the pulse (the origin 
of the zy-coordmate system) to the center of the ellipse, 
and the products at and bt are the semi-minor and semi- 
major axes of the ellipse, respectively. Consequently, the 
parameters a and b are the speeds of expansion of the 
semi-minor and semi-major axes of the ellipse, respec- 
tively, and n is the speed of the center of the ellipse along 
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the x-axis. When rj < a, the ellipse will be incompletely 
"dragged" , and the origin of the pulse will remain inside 
the ellipse at all times. This is the case of a subluminal 
motion of the medium, when the speed of the medium 
u is smaller than c/n. On the other hand, when 77 > a, 
the "dragging" is overwhelming, and the ellipse no longer 
encloses the origin. In this case the medium is moving 
at "superluminal" speeds (c/n < u < c), causing an ex- 
istence of a Mach cone, outside of which no light signal 
can ever reach, and inside of which light is doubled up, 
with two pulses of light in each direction [2^, [32|, [33[ . 

We may translate Eq. into the polar coordinates by 
making the substitutions x — vtcosip and y = vising: 



1 j sin ip 



—5-^ cos If v 

n c 1 



-2U COS If I 1 \r) V + 1l 2 — -^r 

n l n z 



0. 



(7) 



Here, v is the speed of the photon along a light ray di- 
rected at an angle (p from the velocity of the medium 
(see Fig. 1). Equation ([7]) is a quadratic in v, and it can 
be solved for v in terms of the angle tp, the speed of the 
medium it, and the refractive index n of the medium at 
rest. By rejecting the second solution of Eq. © from 
the requirement that v = c/n when u = 0, we obtain: 



v(<p) 



u(l - l/n 2 ) cos ip + [c 2 /n 2 - u 2 (sin 2 ip 



+ cosVn 2 )] 1/2 [l- M 2 /c 2 ] 1/2 } 



1 - u jc (sin ip + cos tp/n ) 



(8) 



From Eq. (j8]) we conclude that the ray velocity in the 
moving medium depends on the angle ip between the di- 
rection of the ray and the velocity of the medium. In this 
sense, an optical medium that is optically isotropic in its 
rest frame will possess an optical anisotropy in the frame 
in which it is moving at a constant velocity. This induced 
optical anisotropy is of a purely relativistic origin, and it 
is different from the usual anisotropy in the crystals [13] • 
The observer in the ^-frame may track the path of a 
light ray in the moving medium by considering it to be 
a "stationary" medium that possesses a relativistically 
induced optical anisotropy characterized by an effective 
refractive index n e f(<p): 



v{ip)' 



(9) 



Recalling Eq. (JSJ) for v(ip), the expression for n e f in Eq. 
© can be re-written in terms of u, n and tp: 



n ef (cp) 



(u/c)(l — n 2 ) cos<p + n[l — u 2 jc 



21 1/2 



2„,2 I T\ V 2 



x [{u/c) 2 (n 2 - 1) cos 2 ip + 1 - n 2 u 2 /c 



x (l-nV/c 2 )" ] 




FIG. 2: The least-time derivation of the refraction formula 
when the material half-space is moving at a constant speed u 
in the direction of the positive a;-axis. 



(10) 



When n = l, Eqs. © and (fTU|) reduce to n e f = 1 and 
v = c for all (p, which describes an expansion of the pulse 
in a vacuum. Also, if u — 0, then n e f = n and v — c/n for 
all tp, which is the case when the medium is stationary. 

If the photon moves parallel to the velocity of the 
medium, n e f k, n±u/c(l — n 2 ) and v w c/n±u(l — l/n 2 ) 
to the first order in u/c, where the sign is taken "+" or 
"-" depending on whether ip = 0° or ip = 180°, respec- 
tively. This longitudinal light drag has been theoretically 
predicted by Fresnel (1818) on the basis of his ether the- 
ory, and experimentally confirmed by Fizeau (1851) in 
his celebrated running- water experiment [32| . 



III. REFRACTION FROM A MATERIAL 
HALF-SPACE MOVING PARALLEL TO THE 
INTERFACE 

A photon traveling in vacuum is incident upon the sur- 
face of a semi-infinite optical material which occupies the 
positive region of the y-axis (see Fig. 2). The material 
half-space is moving at a constant speed u in the positive 
direction of the x-axis, parallel to the vacuum-material 
interface. The material in its rest frame has a refractive 
index n, and is assumed to be homogeneous, isotropic 
and non-dispersive. We will use Fermat's principle to 
derive the formula for the law of refraction of the photon 
from the interface between the vacuum and the moving 
medium. 

Let A be the space-point in the vacuum half-space, 
belonging to the incident ray of the photon, and B is 
the space-point in the moving medium that lies on the 
refracted ray. By (La and ds we denote the shortest dis- 
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tances between the interface and the points A and B, 
respectively, and I is the distance between the orthogo- 
nal projections of the points A and B on the interface. 
According to Fermat's principle, the true path between 
the two fixed points A and B taken by the photon is the 
one that is traversed in the least time. Of all the points 
along the interface, the photon will be refracted at the 
point O whose x-position minimizes the total time tAB 
required for the photon to cover the path from A to B: 



tAB = tAO + toB- 



Here, 



tAO 



(x 2 + dl) 1/2 



(11) 



(12) 



is the time for the photon to travel from A to O at a 
speed of light c, and 



toB 



[(£ - xf + d\ 
v((p) 



(13) 



is the corresponding time from O to B. By v(ip) we 
denote the speed of the photon along the refracted ray at 
an angle tp from the velocity of the medium. If we take 
into account that v{<p) — c/n e f(<p), and use Eqs. (fl~2"]) 
and (I13p . we may re- write the total transit time tAB in 
Eq. dUD as: 



tAB 



(^4) 1/2 + vM [(| _ xf + 4f / 2 



(14) 



We substitute n e f from Eq. (fTT)|) into Eq. fl"3|). use the 
relation 



cos ip = 



l- 



from Fig. 2, and simplify the result to obtain 
(x 2 +d 2 A ) 1/2 (u/c 2 )(l-n 2 ) 



(15) 



tAB 



c 1 — n 2 u 2 /c 2 

„2\V2 ( 



(l-x) 



(n/c) (1 - u 2 /c 2 
1 — n 2 u 2 /c 2 



{u/c) 2 {n 2 - l)(£-xY 



1/2 



+ (l-n 2 u 2 /c 2 )[(e-x) 2 +d 2 B 



(16) 



The principle of least time requires that dtAB/dx = 0, 
and thus: 



= 



x (u/c)(l-n 2 ) n(l - u 2 /c 2 f/ 2 

(x 2 + d 2 ^) 1 ^ 2 l — n 2 u 2 /c 2 l — n 2 u 2 /c 2 

j£~x) J \ (u/c) 2 (n 2 -l){l-x) 2 



d 2 



-1/2 



11 2 2/2 

+ 1 — n u /c 



(17) 



From Fig. 2 we notice the expressions for the angle of 
incidence a and the angle of refraction /3 of the photon: 



sin a 



sin/? = 



{x 2 + d 2 A ) 1/21 

l-x 



[{e-x) 2 + 4F 2 ' 

which we use into Eq. (jTTJ) to get: 

(u/c)(l -n 2 ) _ n(\-u 2 /c 2 f' 2 
1 — n 2 u 2 /c 2 1 — n 2 u 2 /c 2 

sin (3 



(18) 
(19) 



sm a 



[(u/c) 2 (n 2 — 1) sin 2 j3 + 1 — n 2 u 2 /c 2 



1 1/2 



(20) 



If we employ the identity sin/3 = tan/3(l + tan 2 /?) -1 / 2 
into Eq. (|20[) and do some algebra, we obtain a quadratic 
equation in tan/?: 



n ( 1 sin a 



1 - 



tan 2 /? - 



1 



1 — u 2 /c 2 



sin a H — On — 1) 

c 



= 0. 



(21) 



Equation (|21j) has two solutions in tan/3, and the solu- 
tion: 



tan/3 = [(n 2 - l)(u/c) + (1 - n 2 u 2 /c 2 ) sina] 
x < n 2 [l — (m/c) sin a] 2 — (sin a — u/c)" 



x (l-u 2 /c 2 ) 



-1/2 



-1/2 



(22) 



is the only one that is physically correct, which we choose 
from the requirement that when it = it should reduce 
to the usual Snell's law of refraction. 

Equation (]22[) is the formula for the law of refraction 
of the photon from the vacuum-material interface when 
the optical material is moving uniformly parallel to its 
surface. It gives the angle of refraction /3 of the photon 
as a function of the incident angle a, the speed u of the 
medium, and the refractive index n of the medium in its 
rest frame. [Two alternative derivations leading to the 
same refraction formula are given in Ref. [32|.] 

The functional dependence of the angle of refraction j3 
on the incident angle a is plotted in Fig. 3 for n = 1.5 and 
for several different values of the speed u of the medium. 
We immediately notice that the law of refraction in Eq. 
(|2"2")) differs considerably from the usual Snell's law in the 
case of a stationary medium. First of all, when u/c > 0, 
the angle of refraction /3 attain non-zero values even when 
the incident angle a is zero. In other words, the pho- 
ton that enters the moving medium perpendicularly to 
its surface will be dragged by the medium in the mov- 
ing direction. This is the transverse Fresnel-Fizeau light 
drag, theoretically predicted by Fresnel in 1818 and ex- 
perimentally detected by Jones in the early 1970s with 
an accuracy up to the first order of u/c [H, [H, [36l.[37|. 
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FIG. 3: The angle of refraction versus the incident angle 
a for n = 1.5 and for different speeds u of the material half- 
space, when the material half-space is moving parallel to the 
interface. 



The influence of the dragging effect on the refraction is 
also visible for non-zero incident angles, especially when 
they take negative values [38(. In the latter case, if 
u/c > 0, starting from a = 0° down to a certain neg- 
ative incident angle, the angle of refraction /3 is positive 
(see Fig. 3), which means that both the incident and 
the refracted ray of the photon will lie on the same side 
with respect to the interface normal. Moreover, when the 
speed u of the medium becomes greater than c/n 2 , the 
angle of refraction will be positive for any value of the 
incident angle a. 

Equation (J^U) can also be used to calculate the angle 
of refraction of the photon when the medium moves in 
the opposite direction to the one in Fig. 2, that is, in the 
negative direction of the £-axis. However, the situation 
when the photon is incident at a positive angle a upon the 
medium which moves at a speed — u is physically identical 
to the situation which corresponds to a negative incident 
angle —a and a speed u. 

At the end, we emphasize that the absolute value of 
the angle of refraction j3 in Eq. ([22)) cannot become equal 
or greater than 90°, no matter what the values of a, n 
and u are. In other words, it is not possible to observe 
the situation in which the photon would slide along the 
surface of the medium, or to be reflected back to the 
vacuum half-space. To prove this property, we will try to 
find such a value of a which will turn the denominator 
in the right-hand side of Eq. (|2"2")1 into zero. This is the 
condition for the angle to equal ±90°. Hence, a should 
be such that: 



( u ■ \ 2 ( ■ u \' 
n ^1 -smaj — ysmoe J 



= 0, 



(23) 



which is a quadratic equation in sin a, and has two solu- 
tions: 



(sina)i,2 



n ± u/c 



(24) 



nu/c ± 1 

Having in mind that n > 1 and u < c, we have |(n ± 




FIG. 4: The least-time derivation of the refraction formula 
when the material half-space is moving at a constant speed u 
in the positive direction of the t/-axis. 



u/c)(nu/c ± 1) | > 1, which means that there are no 
real values for the angle a to satisfy Eq. (f23l) . 



IV. REFRACTION FROM A MATERIAL 
HALF-SPACE MOVING PERPENDICULAR TO 
THE INTERFACE 

In the following, we consider the case in which the ma- 
terial half-space in Sec. Ill is moving at a constant speed 
u in the positive direction of the y-axis. The geometry 
of the problem is shown in Fig. 4. Observe that unlike 
the case in Sec. Ill, the vacuum-material interface also 
is in uniform motion, having the same constant speed u 
as the medium in the positive y-direction. 

A photon is emitted from point A in the vacuum half- 
space. At this instant, the vacuum-material interface co- 
incides with the x-axis, and cLa is the shortest distance 
between the interface and the point A. After being re- 
fracted from the moving interface, the photon reaches the 
fixed point B in the material half-space. We denote by 
cLb the shortest distance between the interface and the 
point B at time when the photon was emitted from A, 
and by I the distance between the orthogonal projections 
of the points A and B on the interface. 

Due to the motion of the vacuum-material interface, 
the hypothetical points from which the photon may be 
refracted at the interface will not lie on the same hori- 
zontal. According to Fermat's principle, the photon will 
be refracted from a point O that makes AOB the path of 
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the shortest transit time (see Fig. 4). We denote by tAO 
the time required for the photon to go from A to O, and 
by toB the time from O to B. Observe that when the 
photon reaches point O, the interface has already crossed 
the distance utAo m the positive y-direction from its po- 
sition y = when the photon was emitted from A. 

The total time tAB for the photon to travel the path 
AOB is: 

tAB = tAO + toB- (25) 

From Fig. 4 we have: 

{ct AO f = (d A + ut AO ) 2 + x 2 , (26) 
HP)toB? = (d B - ut AO f + (£ - x) 2 , (27) 

where v((3) is the speed of the photon in the moving 
medium in the direction determined by the angle (3. The 
speed v((3) of the photon follows from Eq. ((8]) if we simply 
substitute [3 for Lp. In order to simplify the derivation of 
the refraction formula, we re- write Eqs. (|26[) and (|27|) as: 

t AO = \ [(d A + ut AO ) 2 + x 2 ] lf \ (28) 

t OB = '^^[{dB-utAof + il-xf] 112 , (29) 

with n e /(/3) being the effective refractive index of the 
moving half-space in this case: 



(u/c)(l - n ) cos [3 + n[l - u/c 



2 ,211/2 



[(Wen? 



2„,2 / 21 1 / 2 



1) cos^ 0+1- n z u z /c 



x (1-nV/c 2 )" 1 



Thus, the total transit time tAB in Eq. (|25[) can be re- 
written as: 



1 



[(d A + ut AO ) 2 +X 2 } 1 ' 2 
Uef{P) [{d B -ut A o) 2 + (t-x) 2 } 1/2 • (•'!!) 



tAB 



If we substitute Eq. ([30]) into Eq. (|31j) , and take into 
account the expression: 



cos (3 



utAO 



\{d B - ut AO ) 2 + {£ ~ x)' 



a/2- 



(32) 



from Fig. 4, we obtain an expression for the transit time 
tAB of the photon as a function of the position x of the 
point of refraction O at the moving interface: 



tAB 



[(d A + ut A o) 



2 1 x 2 ] 1 ' 2 u{l-n 2 ){d B -ut A o) 



+ 



1 — n 2 u 2 /c 2 



We now invoke the principle of least time that the value of 
x should be such that the transit time tAB of the photon 
is minimum. We take the derivative of tAB in Eq. (|33p 
with respect to x and set the result to zero to obtain: 







/ dt A o \ 
V dx ) 



(u 2 /c)(n 2 



1) 



1 — n 2 u 2 /c 2 



fdtAo\ _ n(l- u 2 /c 2 ) 
\ dx J 1 — n 2 u 2 /c 2 



( dtAO \ „ 



1/2 



u2 I 2 

— in 



sin/3 

1/2 



1) cos 2 (3 + 1 



where we employed the relations: 

x 



sin a 



cos a 



x 



2 + {d A + ut A oY 

d A + UtAO 



(d, 



sin/3 = 



UtAO) 

l-x 



I 1/2 ' 



[(£ - x) 2 + (d B - ut AO ) 



2l l/2> 



(34) 

(35) 
(36) 
(37) 



from Fig. 4, together with Eq. (|3^|) and the fact that 
tAO = tAo{x)- To find dtAo/dx, we take the derivative 
of Eq. |26|) with respect to x and obtain: 



2c 2 t A o = 2 ( dA + ut Ao)u 



dt A c 
dx 



(30) 

from which we get: 



/ dt A p \ 
V dx J 



x/(ct A o) 



C [1 - (u/c)(dA + utA0)/{ctA0)) 



2x, (38) 



(39) 



By recognizing from Fig. 4 that sin a — xj(c±Ao) an d 
cos a = (dA +utAo)/{ctAo), we recast Eq. into the 
form: 



( dt A p \ 
\ dx J 



sma 



u cos a 



(40) 



We substitute Eq. gOJ) into Eq. {34]), apply the 
trigonometric identities sin j3 = tan/3(l + tan 2 (3)- 1 / 2 and 
cos/3 = (1 + tan 2 /3) -1 / 2 , and do some algebra to obtain 
a quadratic equation in tan /3: 



n (1 — n 2 u /c ) 
l-u 2 /c 2 



. 2 

sin a 



[1 — (u/c) cos a] 



tan 2 /3 



2n 2 (w/c)sina „ (1 — u 2 /c 2 ) sin 2 a „ . , . 

tan/3 - r ^T^ 2" = - ( 41 ) 



1 — (u/c) cos a 



[1 — (u/c) cos a] 



Equation (|41|) has two solutions in tan/3, and its only 
physically correct solution is 



+ 1 



1 1/2 

J [( d B~ut AO ) 2 + (( -x) 2 ] \ .(33) tan/3 = (1 - u 2 /c 2 )sina x [n 2 (7i/c)[l- (u/c) cos 
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FIG. 5: The angle of refraction (3 versus the incident angle a 
for n — 1.5 and for different speeds u of the material half-space 
moving in the positive direction of the y-axis. The portions of 
the curves corresponding to a > arccos(tt/c) are drawn with 
dashed lines. At incident angles corresponding to the dashed 
curves, the refraction does not occur. 



[1 - (u/c) cos a} 2 - (1 - u 2 /c 2 ) sin 2 



1/2 



(42) 



which is the formula for the law of refraction of the pho- 
ton in this case. We reject the second solution of Eq. 
(|4ip as a physically incorrect one, by using the same ar- 
gument as in the deduction of Eq. (22) from Eq. (21). [A 
derivation of Eq. (|4"2")l via the relativistic velocity trans- 
formation formulas is provided in Appendix A.] 

The obtained refraction formula works both for posi- 
tive and negative values of the incident angle a, as well 
as for both moving directions of the medium. The right- 
hand side of Eq. (|4"2"|) is an odd function of a, and hence 
the change in the sign of a does not introduce physically 
new situation. The path of the photon which is incident 
upon the surface at a negative angle —a will be sym- 
metric with respect to the surface normal, to the path 
of the photon corresponding to a positive incident angle 
a. Therefore, in the following, we consider only positive 
incident angles (0° < a < 90°). 

When the medium moves in the positive direction of 
the y-axis, u/c > and the right-hand side of Eq. (|4*2"|) is 
always positive, which means that the angle of refraction 
(3 is between 0° and 90° . The dependence of [3 on the in- 
cident angle a for n — 1.5 and for various positive values 
of u/c is shown in Fig. 5. We see that as u/c varies from 
to 1, the curves corresponding to different values of u/c 
are regularly shifted toward the smaller values of (3, which 
is expectable if one takes into account the dragging effect 
of the moving medium. We emphasize that the portions 
of the curves drawn with dashed lines correspond to in- 
cident angles at which the photon will never reach the 
interface. Namely, for a given speed u of the medium in 
the positive y-direction, there exists an interval of values 
for the incident angle a (a max < a < 90°) for which the 
y-projection of the photon's velocity (coos a) is less than 
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FIG. 6: The angle of refraction /3 as a function on the incident 
angle a for n = 1.5 when the material half-space is moving 
perpendicular to the interface, in the negative direction of 
the j/-axis. When — Ki < u/c < the refraction is regular, 
and the angle of refraction /3 of the photon is in the interval 
[0°,90°]. For — K2 < u/c < —Ki, the refraction is: a) regular 
(0 < (3 < 90°), if < a < a c - and b) backward (90 < /3 < 
180), if Qf c < ct < 90°. If u/c < —K2, the refraction is always 
backward. 



or equal to the velocity of the medium. The value of 



follows from: 



cos a„ 



u/c, 



(43) 



and it can be shown that it is exactly the incident an- 
gle at which the curve j3 = (3(a) attains the maximum 
value. Hence, for incident angles belonging to the interval 
[ckmax, 90°], the refraction becomes impossible. Translat- 
ing to the frame in which the medium is at rest, and thus 
taking into account the aberration phenomenon, we have 
a photon "incident" upon the interface at angles larger 
than 90°. 

The situation changes dramatically when u/c < 0. In 
this case, the denominator on the right-hand side of Eq. 
(|42[) can take zero or even negative values, which implies 
that the angle of refraction (3 can become equal or greater 
than 90°. In Fig. 6 we plotted (3 as a function of a for 
n = 1.5 and for various negative values of u/c. What we 
see in Fig. 6 can be summarized as follows. When the 
medium moves in the negative direction of the y-axis up 
to some value of \u/c\ denoted as K\, the refraction seems 
to be regular. In this range of speeds (— Ki < u/c < 0), as 
the absolute value of the speed increases, the curves are 
shifted toward higher values of (3, and (3 remains smaller 
than 90° for each a. When u/c < —K\, a specific value 
for the incident angle a appears for which the angle of 
refraction of the photon becomes equal to 90° . We denote 
this critical incident angle by a c . For a c < a < 90°, 
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the angle of refraction (3 is greater than 90°. By further 
increasing the absolute value of u/c, the critical angle a c 
decreases, and for some value of \u/c\, denoted by K2, a c 
becomes 0°. For u/c < —k-2, the angle of refraction f3 is 
greater than 90° for all values of a. The quantities k\ and 
K2 are functions on the refractive index n of the medium 
in its rest frame, and as n increases, the values of Ki and 
«2 are becoming closer to each other. It can be shown (see 
Appendix B) that n\ = (1+n 2 )^ 1 / 2 and n 2 = 1/n (in our 
case, for n = 1.5, k% « 0.55 and k 2 ~ 0.67). Evidently, 
the value n 2 corresponds to the smallest superluminal 
speed of the medium. 

The existence of angles of refraction greater than 90° 
in the case u/c < implies that the photon, instead of 
moving in the same general direction as before the refrac- 
tion, moves backward. At first sight, it might seem that 
the photon was reflected at the interface instead of being 
refracted, but this is not what had happened (see Fig. 7). 
Once the photon strikes the vacuum-material interface, it 
undergoes a refraction, enters the medium and continues 
to move into the medium at all times. The curious be- 
havior of the refracted photon is due to the fact that the 
photon, while penetrating into the medium, is dragged by 
the overall motion of the medium. This dragging effect 
is a cause for the y-component of the photon's velocity 
in the moving medium to be negative, and hence, the 
refracted photon will move in the same general direction 
as the medium. 

We emphasize that the path of the photon that un- 
dergoes a backward refraction does not coincide with the 
path which it would trace if it were reflected from the 
interface (see Fig. 7). The path of the backwardly re- 
fracted photon is fully determined by Eq. (HH), and the 
path of the reflected photon obeys the Einstein's formula 
for reflection from a uniformly moving mirror (see Ref. 
[28l [39l . l40l ] ) . In the language of modern quantum elec- 
trodynamics, whether the photon will be reflected or re- 
fracted at the inerface is a matter of probability. In this 
sense, if a beam of photons is incident on the moving 
interface at an angle a, and if the conditions for the ap- 
pearance of the backward refraction are fulfilled, some of 
the photons will be reflected, and some backwardly re- 
fracted. Therefore, in a more realistic view in which the 
reflection is considered also, two different beams spread- 
ing backward will be identified. In this way, we have a 
full consistency between the observations with respect to 
both inertial frames, S and S'. Namely, for the observer 
in S"-frame to which the medium is at rest, both pro- 
cesses are developing in an ordinary way: the refracted 
photon moves forward, obeying the Sncll's law of refrac- 
tion, while the reflected one moves backward, according 
to the usual law of reflection: incident angle equals the 
reflected angle. Hence, the observer in S' clearly distin- 
guishes two rays, and consequently, two rays must exist 
with respect to the observer in the S'-frame: one reflected, 
and one backwardly refracted. 



V. CONCLUDING REMARKS 

We have presented a novel approach to the problem of 
refraction of a light ray at an interface between two homo- 
geneous, isotropic and non-dispersive transparent opti- 
cal materials in uniform rectilinear motion. Our method 
amalgamates the original Fermat's principle and the fact 
that an isotropic optical material at rest becomes opti- 
cally anisotropic if it is moving at a constant velocity. 
The derivation is in the framework of basic optics and 
special relativity, requiring some knowledge of calculus 
at an elementary level. As such, it may be regarded as an 
instructive addendum to the typical introductory physics 
courses. 

We have analyzed in details the refraction at a vacuum- 
material interface, when the material half-space is mov- 
ing uniformly parallel or uniformly perpendicular to the 
interface. In both cases, we considered the plane of in- 
cidence to be normal to the interface and parallel to the 
velocity of the medium. The coincidence of the obtained 
refraction formulas with the ones obtained with the other 
methods confirms the validity of Fermat's principle in the 
presence of uniformly moving boundaries and media. 

In addition, when the material half-space moves per- 
pendicular to the interface, in the direction which is op- 
posite to that of the incident light, we observed a "back- 
ward refraction" of light. To our knowledge, this optical 
effect has not been noticed before in the literature. 

We emphasize that the refraction formulas obtained in 
this paper refer to refraction of a photon, and thus, de- 
scribe the refraction of a light ray. We note that in a set 
of papers on the subject of reflection and transmission of 
a plane electromagnetic wave by a moving medium, start- 
ing with the pioneering papers by Tai and Yeh, the refrac- 
tion formulas for different setups were derived by using 
the notions from the classical electromagnetic theory and 
the Lorentz transformation (4l|, H3, S3, [4^, H^, [H, [4^ • Al- 
though these publications also include the cases discussed 
in our Sections III and IV, the resulting refraction formu- 
las do not match the refraction formulas in Eq. (|22[) and 
P2"| in our paper (48[. This mismatch is due to the fact 
that Tai, Yeh and their successors considered the angles 
of incidence, reflection and refraction of the light wave as 
the angles between the corresponding propagation vec- 
tors and the interface normal. However, we have shown 
that an optical isotropic medium in uniform rectilinear 
motion becomes an optical anisotropic medium. There- 
fore, the propagation vector in the moving medium will, 
in general, not coincide with the light ray in the moving 
medium. In conclusion, the refraction formulas by Tai, 
Yeh and their followers do not refer to refraction of a 
light ray. 

The approach in this paper can be extended to include 
more sophisticated situations, specifically, by taking into 
account the dispersion of the material, or considering var- 
ious three-dimensional cases when the plane of incidence 
is not necessary normal to the interface and not parallel 
to the velocity of the medium. We also suggest a rep- 
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FIG. 7: A space-time simulation depicting the difference between the reflection and the backward refraction. The material 
half-space (the shaded region) is moving in the negative j/-direction at a speed 0.7c. The refractive index of the material 
half-space in its rest frame is n — 1.5. The values of x and y coordinates are given in meters, (a) The situation at time 
t — ns. A photon, located at point (0.2,0.2), is incident upon the vacuum-material interface, making an angle of 30° with 
respect to the interface normal in the vacuum half-space. The interface is located at y = 1.6084478; (b) The path of the photon 
and the position of the interface at time t = 2 ns; (c) The situation at t = 3 ns. The photon reaches the interface at point 
(0.649689,0.97884); (d) The situation at t — 4 ns. The ray br corresponds to the path of the backwardly refracted photon, 
making an angle of 119.974° with respect to the interface normal in the material half-space. The speed of the backwardly 
refracted photon is 0.151c, and its y velocity component is —0.07544c. The ray r is the path of the photon if it were reflected 
from the interface. The angle of reflection is —5.41445° with respect to the interface normal in the vacuum half-space. The y 
velocity component of the reflected photon is —0.99554c; (e) The situation at t — 5 ns; (f) The situation at t = 6 ns. 
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etition of the derivations in Sec. Ill and IV when the 
refraction occurs at an interface between two uniformly 
moving non- vacuum optical materials having different re- 
fractive indices in their rest frames of reference. 
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APPENDIX A: A DERIVATION OF EQ. ([42} BY 
USING THE RELATIVISTIC VELOCITY 
TRANSFORMATION FORMULAS 



of the photon in S-frame are: 



The refraction formula in Eq. (142(1 can also be obtained 
by comparing the path of the photon in the S'- frame 
where the medium is stationary, to the corresponding 
path of the photon in the S'-frame where the medium is in 
uniform rectilinear motion. We first describe the reflec- 
tion with respect to S'-frame. Let the vacuum-material 
interface coincide with x'-axis, with y' > being the re- 
gion of the material half-space, and y' < the region of 
the vacuum half-space. Let the photon be incident on the 
interface from the vacuum half-space in a direction mak- 
ing a positive angle a' with the interface normal. The 
photon is refracted in the material half-space at an an- 
gle (3' with respect to the interface normal, obeying the 
Snell's law of refraction: 



nsin /3', 



(Al) 



where n is the refractive index of the medium at rest. 
Taking that the incident photon moves at a speed of light 
in vacuum c, and the refracted photon at a speed c/n, 
the velocity components of the photon in the S"-frame 
are: 



' !J 2 



csma , 
ccos a', 
(c/n) sin /3', 
(c/n) cos /3', 



(A2) 
(A3) 
(A4) 
(A5) 



where index "1" corresponds to the components of the 
incident photon, and index "2" to the components of the 
refracted photon (for simplicity, we use v' x and v' y instead 
of v' x , and v, ' ,). With respect to S'-frame where the ma- 
terial half-space is moving at a constant speed u in the 
positive y-direction, the photon is incident on the moving 
interface at an angle a at the same constant speed c as 
in the S'-frame, but its speed after being refracted from 
the interface is vp which is a function of the angle of re- 
fraction (3. Thus, the corresponding velocity components 



(A6) 

(A8) 
(A9) 

Dividing Eq. (|A8[) with Eq. (|A9|) . and using the relativis- 
tic transformation formulas for the velocity components 
between S and S': 



csma, 
ccos a, 
vp sin /3, 
vp cos (3. 



<(l- M 2 /c 2 )^ 

l + VyU/C 2 

v', + u 



we obtain: 



tan/3 



1 + v'yu/c 2 ' 



(A10) 
(All) 

(A12) 



Combining Eqs. l(A4jl and (lA5j) with Eq. (|AT|) . we ob- 
tain: 



v 'x 2 = (c/n 2 ) sin a', 

v' y2 = (c/nXl-smV/n 2 ) 1 ^ 



which we substitute into Eq. (|A12|) to get: 
(c/n 2 ) sin a'(l -u 2 /^) 1 / 2 



tan (3 



(c/n)(l-sin 2 a'/n 2 ) 1 l 2 + 



Using the velocity transformation formula: 

, Vx (i-u 2 /c 2 y/ 2 



1 — VyU/C 2 



and Eqs. 



since 



6|) and (|A7[) . we get: 

sina(l -u 2 lc 2 fl 2 
1 — (u/c) cos a 



(A13) 
(A14) 



(A15) 



(A16) 



(A17) 



We substitute Eq. (|A17|) into Eq. (|A15[) and simplify 
the result to obtain the law of refraction of the photon 
inEq. (42|). 



APPENDIX B: CONDITIONS FOR THE 
APPEARANCE OF THE BACKWARD 
REFRACTION 

Because a backward refraction occurs only for u/c < 0, 
we rewrite Eq. (|4"2"|) by putting — k instead of u/c (k = 
\u/c\): 

tan/3 = (1 — k 2 ) sin a x I —nn 2 (1 + Kcosa) 



[1 + k cos a] 2 — (1 — k 2 ) sin 2 a 



1/2- 



(Bl) 
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The photon will be refracted at (3 = 90° when a = a c , 
which is the case when the denominator in the right-hand 
side of Eq. (|Blj) equals zero: 



1/2 



a c = 0°, respectively, and can be obtained from the fol- 
lowing equations: 



(1 - k 2 ) [1 - n 2 (l - nln 2 )] 1/2 - nln 2 (l - K 2 n 2 ) , 



n 2 [1 + k cos a c ] — (1 — k 2 ) (1 — cos 2 a c ) 



(B7) 



Kn 2 (l + KCOsa c ). (B2) [l - n 2 (l - K 2 n 2 )] 1 



/a _ 1 - «2 [1 - » 2 (l - «|" 2 )] 

1 - ft 2 



The last equation is an irrational equation in cosa c , and 
its solution is equivalent to the solution of the system: 



1 — k -\- K 2 n 2 (1 - K 2 n 2 

+2ku 2 (l — K 2 n 2 ) cosa c (33) 
+1 - k 2 - n 2 (l - K 2 n 2 ) =0 



which simplifies to: 



(cos a c )i,2 



cos ei r > 



kti (1 + k cos a c ) > 

{-Kn 2 (\- K 2 n 2 )±{l- k 2 ) 

x[l-n 2 (l-ft 2 n 2 )] 1/2 } 

x [l-K 2 + K 2 n 2 (l-/s 2 n 2 )] _1 
-1/k 



(B4) 

The negative values for cosa c imply a c > 90°. Since 
these solutions are physically unacceptable, we reject 
the solution with the "-" sign in the first equation of 
the system (|B4I) . because this solution corresponds to 
cosa c < 0. What is left to be further analyzed is the 
expression: 



cosa r 



j-ftn 2 (l - K 2 n 2 ) + (1 - k 2 ) 
\l-K 2 n 2 )] 1/2 } 



X 1 — 71 (1 — . 



x 1 1 — K + ft 2 n 2 (l 



2 2 

k n 



)] \ (B5) 



The right-hand side in Eq. ()B5[) is a real number if the 
expression under the square-root is non-negative, i.e. 1 — 
n 2 (l - K 2 n 2 ) > 0. Hence re > (n 2 - X) 1 ! 2 jn 2 = K min . 
Specifically, when k = K m i n , Eq. (|B5|) reduces to cosa c = 
-{n 2 - l) 1 / 2 /™ 2 . Since 0° < a c < 90°, the right-hand 
side in Eq. (|B5j) should satisfy the inequality: 



< 



{-Kn 2 (l - K 2 n 2 ) + (1 - K 2 ) 

x[l-n 2 (l- K 2 n 2 )] 1/2 } 

x [1- k 2 + K 2 n 2 {l- K 2 n 2 )]^ < 1. (B6) 



For a fixed n, one can determine the interval of values 
for k, k £ [ki,K2], such that for each k within this in- 
terval, Eq. (|B5|) has a physically acceptable solution for 
a c . This is the critical incident angle at which the back- 
ward refraction actually appears, and for a < a c the 
refraction is regular (forward), while for a > a c the re- 
fraction is backward. The ends of the interval [ki,K2] 
correspond to the critical incident angles a c = 90° and 



(B8) 



The last two equations are irrational equations in k\^- 
To avoid cumbersome calculations, we substitute x' = 
n 2 {l - K 2 ^ 2 ) in Eq. flE} and x" = 1 - n 2 (l - K 2 2 n 2 ) in 
Eq. (|B8j) to obtain: 



{l-x'f 2 KlX 



1-K 2 ' 
„l/2 = 1 - K 2 X" 
I — K2 



(B9) 
(B10) 



Equation (|B9[) is equivalent to the following system 
kIx' 2 + (l - k?) x' - 1 + k\ = 



x' > 



while Eq. (|B10|) corresponds to: 



k\x" 2 



1 + Ki 



< 





l/«2 



(Bll) 



(B12) 



The quadratic equation in Eq. (|B11[) has two solutions in 
x': x[ = 1 — k\ and x' 2 = —(1 — k 2 )/2k 2 . The solution x' 2 
is not consistent with the inequality in Eq. (|B11|) . since 
x' 2 < (we have taken into account that < k < 1). 
Thus, the only solution of the system (|B11[) is x[ = 1 — k 2 , 
leading to the equation n 2 (l — k\t?) — 1— k 2 , which gives: 



Ki = (1 + n 2 ) 



'1/2 



(B13) 



Similarly, by solving the quadratic equation in Eq. (|B12D 
one obtains two solutions in x" : x'{ = 1 and x 2 — 1/k 2 . 
Since the solution x 2 does not satisfy the inequality in 
Eq. (|B12[) . we reject it and consider only the solution 
x'l = 1, which leads to 1 — n 2 (l — K 2 n 2 ) = 1, and hence: 



k 2 = 1/n. 



(B14) 



Equations (|B13|) and (|B14|) determine the length of the 
interval [k\ , k 2 ] and its dependence on n. We can now 
summarize the results. For k < (1 + n 2 ) -1 / 2 , the refrac- 
tion is regular (forward) for any incident angle a. When 
(1 + n 2 ) -1 / 2 < k < 1/n, there exists a critical incident 
angle a c , such that for a in the interval 0° < a < a c the 
refraction is regular (forward), while for a c < a < 90° 
the refraction is backward. The angle a c can be calcu- 
lated from Eq. (|B5[) . It depends on k, and tends to 
zero as k approaches 1/n. When k = 1/n, a c = 0°, 
and the right-hand side of Eq. (|B1[) becomes an unde- 
termined expression of type 0/0, which can be resolved 
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by using L'Hospital's rule. When k > l/n, the denom- 
inator in Eq. (|B 1 1) is negative for any a. On the other 
hand, the numerator in Eq. (jBljl is always non-negative. 
Therefore, for k > l/n, j3 > 90° for each a. Hence, 
when the medium is moving in the negative direction of 



the y-axis, there exists a critical speed for the medium 
\ u d = K 2 C = c/n, corresponding to the treshold for the 
"superluminal" speeds. Beyond this critical speed, the 
photon will always be backwardly refracted. 
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